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In quantum mechanics, the results of physical measurements performed on a 
physical system prepared in a certain state are random. This fact provides a 
large variety of problems and mathematical and applied statistics offer many 
methods to answer them. For example, statisticians are asked to decide which 
measurement to perform in order to collect most informative observations or to 
rebuild the subjacent quantum object from the resulting observations. In partic- 
ular, physicists have created laser in view of many applications (medical, . . . ). 
An important question is to know whether they have generated the desired one. 
In order to answer this question, we propose a testing procedure based on the 
indirect observations they gathered by a measurement called quantum homo- 
dyne tomography (QHT). To describe our framework and to motivate our work, 
we present in Section 1 the needed physical background. In Section 2, we in- 
troduce our statistical model and define the nonparametric class containing the 
density matrices we are dealing with. To make clear to physicists the statistical 
part, we recall some basic notions on the statistical tests in Section 2.2. The 
testing procedure and the theoretical results assessing its asymptotic behavior 
are presented in Section 3. Section 4 contains an illustration of the proposed 
methodology on simulated data sets while the proof are postponed to Section 5. 

*The authors acknowledge the support of the French Agence Nationale de la Recherche 
(ANR), under grant StatQuant (JC07 07205763) "Quantum Statistics". 
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1. Physical background 

Section 1.1 is a short introduction to quantum optics and we introduce in Sec- 
tion 1.2 the very useful and meaningful pattern functions. In Section 1.3, we 
present some examples of quantum states actually created in laboratory, on 
which our testing procedure may be applied. 

1.1. Short introduction to quantum optics 

In quantum mechanics, the quantum state fully describes all aspects of a physical 
system. Mathematically, the quantum state is described via density operators 
p on a complex Hilbert space TL such that 

1. p is self-adjoint (or Hermitian): p = p* , where p* is the adjoint of p, 

2. p is positive: p > 0, or equivalently (ijj,pip) > for all ijj £ TL, 

3. p has trace one: Tr(p) = 1. 

Moreover, to each measurable physical property or quantity corresponds a self- 
adjoint operator, say X, on the space of states TL. This operator is called an 
observable. Unlike in classical mechanics, when performing a certain measure- 
ment of an observable X of a physical system prepared in a quantum state p, 
the result is in general random and is described by a probability distribution of 
a random variable X. 

In this paper, the studied quantum system is a monochromatic light in a 
cavity described by the state of a quantum harmonic oscillator. Hereafter, we 
consider the complex Hilbert space TL = La(ffi), the space of square integrable 
complex-valued functions on the real line. In quantum optics, the quantum 
harmonic oscillator can be represented by two operators on TL = L 2 (R): the 
position operator Q (or electric field) and the quantity of movement operator 
P (or magnetic field) 



with ipi, ip2 functions of TL. The position Q and momentum P do not commute 
and satisfy the canonical commutation relation: 



where 1 stands for the identity operator on TL. This relation may be understood 
as follows. The measurement of the position Q necessarily disturbs the parti- 
cle's momentum P, and vice versa. They cannot be measured simultaneously. 
Thus, one cannot obtain a couple of random variables (Q, P). However, a linear 
combination 



. dip 2 (x) 
dx 



(1) 



(2) 



[Q,P] =QP-PQ = il, 



X^ = cos(0)Q + sin(0)P, for every phases <j> € [0, n], 
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can theoretically be measured by a technique put in practice for the first time 
in [29] and called quantum homodyne tomography (QHT 1 ), where the phase $ 
is randomly and uniformly distributed on [0,7r]. 

In the case of noiseless measurement and for a phase 4> = 4>, the result would 
be a random variable X with probability density noted p p (-\<p). Moreover, the 
characteristic function of X is 

^lM-MK*) = Tr (pe ttX *) = Tr ( pe «(co S (0)Q+si„(0)P)) _ 

Furthermore, to each state p corresponds a Wigner function W p , which gives 
an equivalent representation of the quantum state p. The associated Wigner 
function W p can be defined rigorously by its Fourier transform ^[Hp] with 
respect to both variables 

W p {u,v) := T 2 [W p ]{u,v) =Tr( ( oexp(mQ + wP)). 

The Wigner function is a mapping from R 2 to R such that // W p (q,p)dqdp = 1. 
Note that the Wigner function can take negative values. For this reason, W p is 
a real- valued function regarded as a generalized joint probability density (quasi- 
probability density) of the two random variables Q and P that we would get 
if we could measure simultaneously the two observables Q and P, which are 
respectively given by equations (1) and (2). It is well-known that its Radon 
transform 3?[W P ] is such that 

p p (x\<f>)=R[W p ](x,<l>). (3) 

Here, 3? denotes the Radon transform, taking functions W p (q,p) on R 2 into 
functions 3?[W p ](a;, </>) on R x [0, tt] formed by integration along lines of direction 
4> and distance x from the origin, expressed as 

/oo 
W p (x cos — t sin (f), x sin <p + t cos 4>)dt. 
-oo 

As 3> is chosen uniformly and independently on [0, n], we can define the proba- 
bility density function p p (x, <f)) of (X, $) w.r.t. the measure —A, where A stands 
for the Lebesgue measure on R x [0, tt] by 

p p {x,4>) = -n[w p ]{x,4>). (4) 

1.2. Fock basis and pattern functions 

An important equation in physics, especially in quantum mechanics, is the 
Schrodingcr equation. The time-independent Schrodinger equation is written 
as 

V> = wV, w G R, (5) 



1 d 2 

~ 2 dx 2 + 2 



1 We refer the interested reader to the book [18] or the paper [16] for further details on 
QHT. 
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where w is the energy level. It turns out that the energies are "quantized" , and 
may only take discrete values noted Wk = k + 1/2, k € N. For a given frequency 
Wk, there are two fundamental solutions: ipk and ipk- One is a normalized func- 
tion, the function ipk, which is such that J ip\ = 1 and is called the regular wave 
function. The other one, the function ifk, is called the irregular one as it cannot 
be normalizable as ipk is. 

The functions {ipk}k^n form a orthonormal basis of L2(R). This particular 
basis, physically very meaningful, is called the Fock basis and is written as 
follows: 

Mx) ■= , * H k (x)e-*'' 2 . (6) 

Here, Hk(x) := (— l) k e x |^e _I denotes the Hermite polynomial of degree k. 

The density operator defined in the previous section corresponds to a density 
matrix under some orthonormal basis. In the Fock basis, the matrix entries p^ k 
of the state p can be expressed as expected values of functions Fj,k(X£,§i) = 
/j,fc(^)e~ l ( fc ~- 7 )*- { , where fj^ — fk,j arc bounded real functions called pattern 
functions [18]. For all j, k £ N, 

Pj,k = JJ o p P {x,ct>)f^{x)e-^ k -^d<t>dx, 

where p p (x, </>) is the joint probability density of {X, $). In other words 

Pjtk = E p [F jtk (X,*)]. (7) 

Equation (7) expresses the idea that one can reconstruct any density matrix 
element pj : k using the pattern functions. First introduced in [20], the pattern 
functions fj t k are well known in physics and are defined in [19] as the first deriva- 
tives of products of the two fundamental solutions ipk and tfk of the Schrodinger 
equation given in (5) for j > k 

d 

fj,k( x ) = WjWVkW) ■ 

A concrete expression for their Fourier transform using generalized Laguerre 
polynomials can be found in [28, 4]. For j > k 



f jtk (t) =n(-i)i-^ 2 ±^\ t \V- k e- t lLi-\^), (8) 

where fj t k denotes the Fourier transform of the pattern function fj tk and L^(x) 
denotes the generalized Laguerre polynomial. We note that the pattern functions 
fj,k(x) are even functions for even differences j — k and odd functions for odd 
ones 

/i,*(-aO = i-iy-^^ix). 
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Table 1 
Examples of quantum states 

Vacuum state 

• P0,0 = 1 res t zero, 

• Pp(x\4>) = e~ x /0r. 

Single photon state 

• pi.i = 1 rest zero, 

» p P (x\4>) = 2x 2 e~ x /yjr. 

Coherent-qo state go G K 

• Pj,k = e-\"o/^\ 2 (q /y/2y+ k /JjW., 

• p P {x\<t>) = exp(-Q - g cos(</>)) )/y / n- 

Squeezed state (A/, 8) G R^_, £ G K 

. p,, fc = C(M,0(5tanli(0)' s+;/ ^(5)fffe(<5)/Vi!fe!. 

• p p (x|0) = exp [(sin 2 (0)(xe -2 £ cos(c/>) — (xcos(^)) — a)e 2 £) 2 ) / (e 2 ^ sin 2 (0) + e~ 2 ^ cos 2 (<ji)) 

X — e~ 2 £(a'cos(</>) — a) 2 — e 2 ^x 2 sin 2 (</>)j / ^(e 2 ^ sin 2 (^) + e~ 2 £ cos 2 ((j>)). 
Thermal state f3 > 

• Pj.fc = (1 " e^V^l^, 

» p P Qz#) = x /tanh(/3/2)/7rcxp(-x 2 tanh( / 3/2)). 

Schrddinger cat-go go > 

• pj.fc = 2(<jo/v / 2) :,+fe / ^-\/i!fc!(exp(ijQ/2) + exp(— gg/2))J , for j and fc even, rest zero, 

• pp(x\<f>) = (exp(-(:r - go cos(^)) 2 ) + exp(-(x + go cos(<j>)) 2 ) 

+2cos(2g xsin(^))exp(-x 2 - g 2 cos 2 Q))) / (2^(1 + cxp(-g 2 ))) . 

1.3. Examples of quantum states 

We present in Table 1 examples of pure quantum states, which can be created at 
this moment in laboratory and belong to the class lZ(B,r, L) with r = 2. A state 
is called pure if it cannot be represented as a mixture (convex combination) of 
other states, i.e., if it is an extreme point of the convex set of states. This is 
equivalent to the density matrix being a one dimensional projector, i.e., of the 
form p = P^f, for some unit vector ip. Equivalently, a state p is pure if Tr(p 2 ) = 1. 
All other states are called mixed states. 

Let us discuss these few examples of quantum states. Among the pure states 
we consider the single photon state and the vacuum state, which is the pure state 
of zero photons. Note that the vacuum state would provide a random variable 
of Gaussian probability density p p {x\4i) via the ideal measurement of quantum 
homodyne tomography. We consider also the coherent-q$ state, which charac- 
terizes the laser pulse with the number of photons Poisson distributed with an 
average of M photons. The Squeezed states (see e.g. [7]) have Gaussian Wigner 
functions whose variances in the two directions have a fixed product. The param- 
eters M and £ are such that M > sinh 2 (£), C(M, £) is a normalization constant, 
a = ((M - sinh 2 (£)) 1/2 )/(cosh(£) - sinh(£)), and 6 = (a/(sinh(20)) 1/2 . The 
Schrddinger cat state is described by a linear superposition of two coherent vec- 
tors (see e.g. [25]). Table 1 gives some explicit density matrix coefficients 
and probability densities p p {x\<j>). 
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2. Problem formulation 



In this paper, the studied quantum system is a monochromatic light in a cavity, 
whose state is described by an infinite density matrix p on the Hilbcrt space 
TL = L2(M). In this setting, a convenient representation of a quantum state can 
be obtained by the projection onto the orthonormal Fock basis ("0fc)fc6N defined 
in (6). In quantum optics, physicists produce quantum state of light and via 
QHT measurements, they gather independent identically distributed random 
variables containing information on the unknown, underlying quantum state p. 
In an ideal framework, the results of measurements would be (Xg, $i)t=\,...,ni 
independent identically distributed random variables with values in K x [0,7r], 
with p p (x, 4>) the probability density function of (X\, $i) defined in(4). 



2.1. Statistical model 



In this paper we consider a more realistic model in presence of an additional 
independent Gaussian noise. In practice from n identical, independent QHT 
measurements, we do not collect data (Xe,&e), but we observe (Ye, &i)i=i,..., n 
independent identically distributed random variables, as error is add such that 



Y t := y/rjXt + ^(1 - T7)/2 fc, (9) 

where £f is a sequence of independent identically distributed standard Gaussians, 
independent of all (Xi, $£). The detection efficiency parameter r], < r) < 1, is 
a known parameter and 1 — r\ represents the proportion of lost photons due to 
various losses in the measurement process. 

First recall that for any functions /, g : R — ► R, we denote by f * g the 
convolution product 

/ * g(y) = f f(y- t)g(t)dt. 



Due to (9), the density p 7 ^ of (Yi, &e ) is given by the convolution of the density 
Pp{'/ ^/Vi4>)l y/V with N n the centered Gaussian density of variance (1 — rj)/2. 
In other terms 



p%y, 0) = ( -^ Pp [ —,4> ) * ) (y), Vy G R, G [0, tt]. 



1 

This corresponds to a severely ill-posed inverse problem as the additive noise is 
super-smooth Gaussian. In the Fourier domain, this relation becomes 

FxWli-Mit) =Fi\p P (;<t>WVv)N r >(t), (10) 

where T\ [p v p {-, <p)\ denotes the Fourier transform with respect to the first variable 
and N v the Fourier transform of the Gaussian density N n . 

We suppose that the unknown state belongs to a natural class of states from 
the point of view of applications, the class TZ(B,r,L) for B > and r G]0, 2], 
defined by 

K{B, r, L) := {p quantum state : \p jjh \ < Lexp(-B(j + k) r/2 )}. (11) 
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Let us note that all the states described in Table 1 belong to this class with 
r = 2. 

2.2. Statistical tests 

An important problem in quantum optics is to check whether the produced light 
pulse is in the desired known quantum state r or not. The purpose of this paper 
is to answer this question via goodness-of-fit testing. 

More precisely, we consider here the problem of nonparametric goodness-of 
fit testing from the data (Ye, <&e) for I = 1, . . . , n; i.e. given r £ 1Z(B, r, L), we 
define the null hypothesis H and the alternative hypothesis Hi as follows: 

f H : p = r, 

\Hi(C,<p n ) : p€ll(B,r,L) s.t. that \\p - r|j 2 > C ■ <p n , 

where (p n is a sequence, which tends to when n — > oo. The physical interpre- 
tation of such a test is to check whether the produced light pulse is in a known 
quantum state r (H is accepted), or not (Hi is accepted). Here, the distance 
between the unknown state p and the presumed state r is measured by the 
squared-L2-distance: 

\\P-T\\l= £ K fc -r„ fc | 2 . (12) 

j,k>0 

In nonparametric statistics, different tools are developed to evaluate the ac- 
curacy of a testing procedure. First, let us begin by reminding of some basic 
definitions. There arc two important errors made in a statistical decision process: 

• First-type error: the test will reject a correct null hypothesis, 

• Second-type error: the test will accept a false null hypothesis. 

Given a test procedure fl n such that Q n = when we accept Hq and Q n = 1 
when we reject Hq and decide Hi, we denote by a — P T [fl n = 1] and (3 = 
(3(p) = P p [£l n = 0] the probabilities to make a first-type error and a second- 
type error respectively under r defined in H and p satisfying Hi(C,ip n ), We 
would like to control the sum of these two probabilities and we do it as described 
in Definition 1. 

Definition 1. For a given < A < 1, a test procedure £l n satisfies the upper 
bound (13) for the testing rate ip n over the smoothness class lZ(B,r, L) if there 
exists a constant C* > such that for all C > C* : 

limsup \P T [n n = 1] + sup P p [n n = 0]} < A, (13) 

where P T denote the probability under p = t defined in Hq. 

In the statistics, we may also compare the power of several testing procedures 
having first-type error less than some fixed a and choose the most powerful 
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procedure. The power of a statistical test is the probability that the test will 
reject a false null hypothesis (that it will not make a second-type error) and we 
denote it by II such that 

II = 1 - f3 = P p [fl n = 1], under p satisfying Fi(C, tp n ). (14) 

Remark that, the probability of a second- type error decreases as the power 
increases. In Section 4, to evaluate the performance of our testing procedure 
O n , we estimate empirically the power of our test. 

2. 3. Outline of results 

The problem of reconstructing the quantum state of a light beam has been 
extensively studied in quantum statistics and physical literature. Methods for 
reconstructing a quantum state arc based cither on density matrix or on Wigncr 
function estimation. The estimation of the density matrix from averages of data 
has been considered in the framework of ideal detection [13, 12, 20, 3] as well as 
in the more general case of an efficiency parameter r\ belonging to the interval 
]l/2, 1] (cf. [10, 14, 11, 27]). The case -q e]0, 1] has been recently treated in [4]. 
The latter paper provides also rates of convergence in L2 loss for an estimator 
of the Wigner function. The problem of pointwise estimation of the Wigner 
function has been previously studied in [16] for ideal data and in [9] for noisy 
data. It should be noted that the results of [16] and most part of the results in 
[9] are asymptotically minimax not only in the rate, but also in the constant. 

In the present work, the goodness-of-fit problem in quantum statistics is 
considered. There is a large literature on nonparametric testing procedures for 
the goodness-of-fit of probability distributions. First of all, let us mention the 
family of test procedures built on certain distances between empirical cumulative 
distribution functions (c.d.f.), such as the Kolmogorov-Smirnov and the Cramer- 
Von Mises test statistics, for which extensive results in terms of asymptotic 
efficiency were established [23]. In order to be more sensitive to low- frequency 
components or narrow bumps for powerful discrimination, test procedures based 
on distances between densities were proposed, such as the Bickel-Rosenblatt 
test [G, 2] for the L2-distance and the test of [1] for the Li-distance. Theoretical 
results on such test statistics naturally stems from their nonparametric function 
estimation counterparts. 

In order to compare nonparametric testing procedures, many approaches were 
proposed, as reviewed in [23, 17]. A common approach is to analyze the power 
against sequences of local alternatives {f n }n>i of the form /„ = fa + (f n g where 
g is the direction defining the sequence of local alternatives, and where (p n — > 
as n tends to infinity. Typically, while nonparametric test statistics achieve non- 
trivial power (i.e. the power of the test is strictly larger than the first type error) 
against directional alternatives /„ = fo + f n g when ip n = Cn^ 1 / 2 , they achieve 
nontrivial power against nondircctional alternatives /„ = / + f n 9n only for <p n 
slower than n~ x l' 2 . In other words, achieving nontrivial power uniformly against 
a large class of alternatives comes at the price of a slower rate than for para- 
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metric testing. The minimax distinguishability framework described in a non- 
asymptotic setting in [5, 15] and in an asymptotic setting in [17], allows to give 
precise statements about this phenomenon by characterizing the discrimination 
rate ip n depending on a smoothness index of the class of alternatives one wish 
to discriminate from the null hypothesis. Sharp minimax results with pointwise 
and sup-norm distances have been established in [21] for the regression model 
and in [2(>] for the Gaussian white noise model for supcrsmooth functions. In 
[ ], goodness-of-fit testing in the convolution model have been considered and 
minimax rates for testing in IL,2-norm from indirect observations have been es- 
tablished. The first testing procedure adaptive to smoothness of the alternative 
function was proposed by [30]. 

In quantum statistics framework, for the problem of discriminating between 
two different and fixed states we mention [24] among others. They have establish 
lower bound for the Bayesian error probability. 

The remainder of the article is organized as follows. In Section 3, we present 
our testing procedure and our theoretical results. We study in Section 4 the 
numerical performance of our testing procedure. The proofs are deferred to 
Section 5. 

3. Testing procedure and results 

Our testing method is based on a projection method on pattern functions. We 
describe first the pattern functions adapted to our setting and derive some 
useful properties. We define and study the test statistics and give the main 
results concerning the testing procedure. 

3. 1 . Pattern functions 

To take into account the detection losses described by the overall efficiency n, it 
is necessary to adapt the pattern functions. When n e]l/2, 1], we denote by fV k 
the suitable functions introduced in [4] and defined by their Fourier transform 
as follows: 

f], k {t)-=kk{t)^- (15) 

In this paper we develop a procedure for n G]l/2, 1], but it is possible to get 
quite similar results with the same procedure for < i] < | by using modified 
pattern functions, those introduced in [4]. We restrict our study to the more 
interesting case r\ G]l/2, 1] as in practice r\ is around 0.9. The following lemma 
provides useful upper bounds on the L,2-norm of the fj t k and the Loo-norm of 
the fj k . From now on, we denote 7 := -jj?- 

Lemma 1 (Aubry et al. [4]). There exists a constant C2 such that for N large 
enough 

E ll/i.*lla<^. 

j,k=0 
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For i] e]l/2, 1[ and 7 = -3^, i/iere exists constant swc/i £/ia£ /or iV Zarge 



AT 
J,k=0 



<ClAf-ie 1677V . 



For the Loo-norm of the fj k , this lemma is slightly different from Lemmata 4 
and 5 in [4] where the sum is over j + k = 0, . . . , N. The proof remains similar. 

3.2. Testing procedure 



In this part, we propose a testing procedure which allow to choose among the 
hypothesis Hq and Hi defined in Section 2.2 by 

Ho- P = t, 

#i(C, </?„): p € 1Z(B, r, L) s.t. \\p - r|| 2 > C ■ ip n , 



As in the alternative Hi, the true value of ||p — r|| 2 is unknown, we have to 
derive an estimator of this quantity. Then, in a second time, we provide our 
testing procedure Q n . 

For the known efficiency parameter 77 €]l/2, 1], we define an estimator M n , 
also called a test statistic of \\p — T^.the estimator M n is as an U-statistic of 
order 2 based on indirect observations (Yi, $i)i=i,..., n - For N := N(n) — > 00 as 
n — > 00, M„ is defined by 



1 



N-l n 



n(n — 1) ^— ^ ^— ' 



— ; CP m -T iifc 



(16) 



where Fj k (x,cf) := ' k {x)e^ k ~^ uses the pattern functions defined in (15) 
and a denote the complex conjugate of a. Let us note that the density matri- 
ces are infinite, thus we truncate the sum to N. The parameter N is called 
the bandwidth and has to be optimized. We study the properties of M n in 
Section 3.3. 

Let us discuss the construction of our estimator M n . First remark that we 
do not estimate 



j,k=0 



N-l 
j,k=0 



j,k=N 



IN 



R 



N 



but the truncated part djv- For an appropriate choice of N, the residual term 
Rn should be small due to the considered decrease condition (11). Moreover, 
we have underlined previously that pj_u = E p [Fj t k{X, <£>)], thus we can estimate 
Pj,k by the empirical version p^ k , < j, k < N (see [4]). For the case 77 G]0, 1[ 
we shall use the pattern functions defined in (15). We estimate then djy by an 
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U-statistic of order 2 given by (16). By the plug-in method, d,N can be estimated 

by 

JV-l JV-l 

- E 

j,k=0 



^ n = E \Pj,k - 

3,k=0 

However, it is well-known in statistics that such estimators have a too large bias 
term. 

Now, for a constant C > and some threshold t n > defined later, we can 
define our testing procedure, based on the test statistic M n defined in (16), as 

O n = l(|M n | >C-t 2 n ). (17) 

Here, VL n = when p — r and Q„ = 1 when the distance between the unknown 
state p and the presumed state r is larger than Ct^. We study the upper bounds 
of our procedure Q n in the sense of Definition 13 in Section 3.4. 



3.3. Properties of the estimator M„ 



We first remark that each element of the density matrix pj^k, such that j, k < N, 
is estimated with no bias. Moreover, the estimator M n is unbiased under p = r 
defined in Hq (Remark 1). Remark 2 derives useful tools for the proof of the 
upper bounds in the sense of Definition 1. From now on, we denote by E p the 
expected value under p satisfying Hi(C, ip n ) and E T under p = r defined in Hq. 

Remark 1. Note that forO<k<j<N-l and r) G]l/2, 1] 



F\i — ,* 



Pj t k and E T 



F\,[ — 



Indeed, from Plancherel formula 
E 











-i{k-j)c, 




e -i(.k-j) 


4> 1 t. 


k {t)e^ t2 F^p^-^^Wdt 



From equation (10), and from Plancherel formula 



E„ 



2tt 



f jt k (t)e— 1 T x \ Vp {- , </>)] (i) JV" {t/yfq)d4dt 



-i(k iWf^xjppfafydffxix = p. jM . 



Thus, 



N-l 



E p [M n ] = J2 (Pj,k ~ T j,kf and E T [M n ] = 0. 
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Remark 2. For N large enough and p belonging to the class lZ(B,r, L) 



N-1 

E^ 

j,k=0 



< 



C^n-V^n r/ e]i/2, l[, 
c-c 2 n 17 / & n = i, 



where and Ci are constants defined in Lemma 1 and C is a positive constant. 
Proof. Since £ p $) | 2 ] = E p [\fl k (l=,$)\ 2 ] For „ e]l/2, 1] 



N-1 

E E ? 

j,k=0 



3,k \ yffj 

Y 



P 3,k 



< 



jL^j,k—0 J j,k 



N-1 



j,k=0 



f" 



r?e]i/2,i[, 
v = i. 



For rj = 1, we first apply Lemma 6 in [4] where it has been established that 



p p {x,4>)d(j} < C, 



sup 

xGRJq 



where C is a positive constant. Hence, the result is a direct consequence of 
Lemma 1. □ 

In order to choose the optimal the bandwidth N, we do the classical bias/variance 
trade-off of the risk E[\M n - \\p - t\\1\ 2 ] = B 2 (M n ) + V(M n ), where the bias 
and the variance term are respectively defined as follows: 

B{M n ):= \E[M n ]-\\p-T\\l\, 



V{M n 



E 



\M n - E[M n ] 



From now on, we denote by £> p , V p the variance and the bias terms under p 
satisfying H\{C, (p n ) and V T , B T the variance and the bias terms under p = r 
defined in Hq. Note that under p — r defined in Hq, the bias term is written 
B T (M n ) := \E T [M n ]\ and is equal to 0. Hence, the following Propositions 1 
and 2 evaluate these quantities. 

Proposition 1. For r €]0, 2], and ij G]l/2, 1] we have 
B T (M n ) = 0, 

B p {M n ) < C B N 2 - r ^e- 2BNr/ \ for N large enough , 

where Cb is a positive constant depending only on B, r and L. 

Proposition 2. Forr g]0, 2], rj = 1 and for N large enough such that N 17 ^ 6 /n — ► 
as n — > oo, we have 



V p {M n ) < 



8C-C 2 
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For r s]0,2], r\ e]l/2,l[, 7 = (1 — 7])/(An) and for N large enough such that 
N~ 2 /^e 1&lN /n — ► as n — > 00, we have 

onn at— 2/3 

V p (M n ) < ^£L e 16 ^, 

n 

(rV ^2 AT— 4/3 

V T (M n ) < e 32 ^. 

The constants and C2 are defined in Lemma 1 and C is positive constant. 



3-4- Main results 

In the following Theorem, we establish upper bounds for the testing rates in 
the sense of Definition 1 where p is the unknown density matrix supposed to 
belong to the class lZ(B,r,L) defined in (11). Theorem 1-(1) deals with the 
ideal detection case while Theorem l-(2) and Theorem 2 take into account the 
Gaussian noise. 

Theorem 1. The test procedure fi n defined in (17) for the bandwidth N(n), 
the threshold t n and the constant C* satisfies the upper bound (13) for the rate 
tp n such that 

1. for r s]0, 2], r] = 1, the bandwidth N(n) := Ni is equal to 

/bgn (log log n) 2 ^ 

\ AB AB J ' v ' 

and the rate 

ip 2 n =t 2 n = n- 1 /2(i gn)^, (19) 

2. for r = 2, r) s]l/2, 1[ and 7 := -j-^, the bandwidth N(n) := N2 is equal to 

log(n) / 8 1og(logn) \ 
N2 -4(A 1 + B) ( v 1+ 3^oi(^rJ' (20) 



and the rate 



n 1 12-y-B B 

<Pn = hi = l°g(«) 3(4 ^+- B) U 2(4 7 + B) _ (21) 



Theorem 2. For r e]0, 2[, 77 e]l/2,l[ and 7 = ij- 2 , i/ie test procedure fl n 
defined in (17) for the bandwidth N := 7V3 solution of the equation 

I67N3 + ABNl /2 = log ra, (22) 
i/ie threshold t n and the constant C* satisfies the upper bound (13) for the rate 

<£n 

^=4 = iV 3 2 - r /V 2 *< 2 . (23) 
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We can remark that 



limsup \ P T [n n = 1] + sup p p [n n = 0] \ = 




Theorem 1 and 2 provide upper bounds for our testing procedure defined in 
equation (17) for the testing rates <p n . We first remark that our procedure gives 
nearly parametric rate up to a logarithmic factor when we are in the framework 
of ideal detection (Theorem l-(2): no noise) and supersmooth corresponding 
Wigner functions for all r e]0, 2]. In the setting of Theorem l-(2) and Theorem 2, 
there are good reasons to believe that our testing procedure achieves optimal 
rates. This remark follows from a recent work of Butucea ([8]), where the author 
establishes minimax rates for testing in L,2-norm from indirect observations. 
However, we do not attempt to go that far in this paper. 

In the framework of quantum statistics, one can investigate another testing 
approach based on kernel type estimator of the Wigner function describing 
the quantum state. Such a testing procedure can been directly derived from 
the kernel estimator in [22] of the quadratic functional JfWp, where W p is 
the Wigner function associated to the quantum state p. Our test problem is 
equivalent to the following: 



where ip n is a sequence which tends to when n — > oo. We conjecture that its 
performances are comparable to those found in this paper and we will leave this 
analysis for a separate work. 

4. Simulations 

From now on, we set r = 2. The purpose of this section is to implement our 
testing procedure and to investigate its numerical performances. Our motivation 
is, given a density matrix r <E 1Z{B,r), to decide whether Hq or Hi is accepted 



where <p n is a sequence, which tends to when n — > oo. 

We propose to simulate two different situations. In the first one, case A, we 
consider quantum states easily distinguishable, while in the second one, case 
B, we deal with quantum states, which are quite similar and it is difficult to 
differentiate between them. For r defined in Hq, we sample our procedure from 
different density matrices p satisfying Hi(C,<p n ) such that in 

• the case A: r is the vacuum state, while 

— a) p is the vacuum state (p = r), 

— b) p is the single photon state, 

— c) p is the Schrddinger cat-3 state, 




'n 



{ 



H : P = t, 
Fi(C,0: p e TZ(B, r, L) s.t. \\p-T\\ 2 >C-<p. 
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Figure 1. Density matrices of the vacuum state, the single photon state, the Schrddinger 
cat-3, the coherent-3 state and the coherent-%/6 state respectively. 

• the case B: r is the coherent-3 state, while 

— a) p is the coherent-3 state (p = r) , 

— b) p is the coherent-\[% state, 

— c) p is the Schrddinger cat-3 state. 

The latter case is a more complicated situation for two main reasons. The 
Schrddinger cat-qo state corresponds to a linear superposition of two coherent 
states (±(7o)- Moreover, the probability density of a coherent-qo state is Gaussian 
with a mean proportional to qo (see Table 1). Figure 1 represents the density 
matrices of the states we consider in our simulations. In each situation (cases 
A and B), we shall distinguish the case r\ = 1, which corresponds to the ideal 
detection, from the case r\ = 0.9, i.e. when noise is present. The latter case is 
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the practical one in laboratory. From now on and for both situations 

• when r/ = 1. we set N = 15, 

• when r} = 0.9, we set N = 14 and N = 13. 

The values of N have been chosen after different simulations for different values 
of N. We notice that, when we are in presence of noise (77 = 0.9), we have to 
take N smaller than the one we choose in the ideal setting. It can be explained 
by the fact that the variance term of our estimator increase exponentially with 
N in the case 77 = 0.9 (see Proposition 2). 

To implement our procedure, we first compute our modified pattern function 
fj k (x) in Section 4.1. We implement in Section 4.2 the estimator M n defined in 
(16) and finally study the performance of our test procedure f2„ defined in (17) 
in Section 4.3. 

4-1- Pattern functions fj k 

To implement our procedure, we need the modified pattern function fj k (x) 
defined in (15) for all j > k. For this purpose, we compute here the inverse 
Fourier Transform fj k (x) of the fj k (t) given explicitly by (15) and the gener- 
alized Laguerre polynomials. Previous authors have used a different method to 
implement the pattern functions in [19, 3] via the following recurrence relation 
given in [18]: 

fj,k( x ) = Zxipk^Pjix) - V 2 ( k + l)^k+i{x)^j{x) - y/2(j + l)ip k (x)ipj+i(x), 

for j > k, otherwise fj,k(x) = fk,j(x)- We display in Figure 2- (a) the corre- 
sponding graphical representations of the pattern functions up to a constant ir 
as in physical literature 7r _1 fj. k instead of fj^k are often called pattern functions. 
Figure 2-(b) represents some modified pattern functions 7r -1 /?fc for i] = 0.9. 




-5 5-5 5 -5 5-5 5 



(a) fj, k = fl k for V = 1 0>) /j; k , 17 = 0.9 

Figure 2. Examples of pattern functions. 
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For some of fj^k we expressed below their explicit form. 
fo,o(x) = 2 - 2e-* 2 aV^Erf i(x), 



hA x ) 



30 



-2e a;2 .T(-3 + 2a; 2 ) + v^r(l - 8a; 2 + 4x 4 )Erf i(x) , 

2e ;E2 (-4 + 27a; 2 - 24a; 4 + 4a; 6 ) 

x 6 ) Erfi(x)] , 
2e x2 (-30 + 435a; 2 - 865a; 4 + 526a; 6 - 116a; 8 + 8a; 10 ) 



2^3 

-V^x (21 - 74a; 2 + 52a; 4 



-y/^x (225 - 1425a; 2 + 2160a; 4 - 1160a; 6 + 240a; 8 - 16a; 10 ) Erf i(x)] . 



4-2. Implementation of M n 

The purpose of this section is to investigate the performance of the estimator 
M n of ||p - -7-111 defined in (16) by 



M n 



1 



JV-l 



n(n - 1) Z-i 2^2 

K ' j,k=0 i=£m=l 



$e - t. 



3,k 



r 3,k 



where F^ k {x,4>) = fjk( x ) e ^ anc ^ T 1S defined in H . Table 2 gives the 
actual values of \\p — r|j 2 in all considered setups. 

From now on and for all considered cases (A and B), the procedure to sim- 
ulate one value of M n is designed as follows. We first sample n = 50 000 i.i.d. 
data from each of the three probability densities p p defined in Table 1 for each 
considered p (Figure 1). The Gaussian noise component £ is simulated indepen- 
dently with a variance equal to ^(1 - rj)/2. Then, we run independently of the 
other runs, 1000 values of M n denoted by {Af*},. 100Q , where each value 

M% is based on n = 50 000 i.i.d. noisy data (Yt, <&i). 

We recall that the choice of TV = 15 corresponds to the setting of ideal 
detection with r\ = 1, while N = 14 and N = 13 deal with the noisy detection 
with r\ = 0.9. Note that for p = t we expect the M„ to be close to 0, while for 
p satisfying H\{C, ip n ) we expect the M% to be close to \\p — r||| > 0. 

In frameworks A and B, we draw boxplots of the {M*},, 1 0Q0 in Figures 3 
to 8. For a clearer reading, these boxplots (Figures 3 to 8) have been moved in 
an appendix at the end of the paper in Section 6. We summarize the median 
values in Table 3 and we can make the following remarks. 



Table 2 

The values of \\p — t||| for r the vacuum state (case A) and r the coherent-?, state (case B) 



Case A; 


a) p vacuum 


b) p single photon 


c) p Schrodinger cat-3 







2 


1.9556 


Case B: 


a) p coherent-3 


b) p coherent-\[§ 


c) p Schrodinger cat-3 







0.2812 


0.9999 
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Table 3 

The median values of Af„ for r the vacuum state ( case A ) and r the coherent-?) state 

(case B) 



Case A: 


a) p vacuum 


b) p single photon 


c) p Schrodinger cat-3 


N = 15, n = 1 


0.0043 


1.9972 


1.9417 


N = 14, n = 0.9 


0.0256 


1.9936 


1.9360 


N = 13, n = 0.9 


0.0144 


1.9943 


1.9378 


Case B: 


a) p coherent-3 


b) p coherent-^/Q 


c) p Schrodinger cat-3 


N = 15, r) = 1 


-6.8721 ■ 10~ 4 


0.2688 


0.8962 


N = 14, ij = 0.9 


-0.0466 


0.1704 


0.7735 


TV = 13, J7 = 0.9 


-0.0376 


0.1847 


0.7840 






Table 4 




77ie empirical values 


of MSE for t the 


vacuum state ( case A ) and r the coherent-3 state 






(case B) 




Case A: 


a) p vacuum 


b) p single photon 


c) p Schrodinger cat-3 


AT = 15, ?7 = 1 


2.3988 ■ 10~ 5 


1.3175 • 10~ 4 


4.7877- 10" 4 


AT = 14, 77 = 0.9 


0.0181 


0.0056 


0.0071 


N =13, n = 0.9 


0.0069 


0.0025 


0.0051 


Case B: 


a) p coherent-3 


b) p coherent-^/Q 


c) p Schrodinger cat-3 


AT = 15, r} = 1 


1.4303 ■ 10~ 6 


1.9473 ■ 10~ 4 


0.0110 


N = 14, ?7 = 0.9 


0.0044 


0.0132 


0.0531 


AT = 13, ?7 = 0.9 


0.0019 


0.0093 


0.0484 



• In the setting A: the procedure M n is a good estimation of \\p — t| 1 1 both 
for r) = 1 and 77 = 0.9. 

• In the setting B: when p — r, the procedure M n shows excellent results 
since the M% are close to \\p— t\\\ = when either rj = 1 or 77 = 0.9. In the 
ideal case and when p is the coherent-\/6 state, the procedure M n gives a 
good result since the median of the boxplot of the M% is equal to 0.2688 
and the true value \\p — t\\^ = 0.2812. Otherwise, the procedure M n under 
evaluates the distance \\p — t|| | . 

In order to evaluate the quality of our procedure M n we estimate the mean 
square error MSE = E p [\M n — \\p— 7-| 1 2 1 2 ] as the average over the 1 000 indepen- 
dent runs of \M n — \\p — tH^ 2 - In other words the MSE is empirically assessed 
as 

1 1 000 

^EWMiP-niil 2 . 

fe=i 

with — t\\\ given by Table 2. Table 4 summarizes the results. As we have 
already noticed the procedure M n gives excellent results in every cases but it 
has a larger MSE when we evaluate the distance between the coherent-3 state 
and the Schrodinger cat-3 state: MSE is equal to 0.0531 and 0.0484. 
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4-3. Studies of the performance of our test procedure Vt n 

In this part, we would like to confirm the performance of our testing procedure 
f2„. In order to appreciate it, we are interested in the power of our test II, such 
that II = -P p [il„ = 1] under p satisfying H\(C,ip n ). In this view, we want to 
compare our estimator M n with a threshold v n (with v n = C*t 2 n in (14)) and 
decide as follows: 

• if \M n \ > v n , we accept Hi, 

• otherwise we accept Hq. 

In the statistical literature, the parameter v n = v n {a) is also called the critical 
value associated to a = P T [il n = 1], which is the probability error of first type 
(under t defined in Hq) and defined in Section 2.2. In our framework, we set 
a = 1% and a = 5%. As we don't know the density probability of M n under 
p = t defined in Hq, we shall evaluate empirically the threshold v n at the 
testing level a for our particular choices of r, n, a and rj. In this purpose and 
independently of the future runs, we first compute, for n = 50 000 and for all 
considered cases 1 000 independent values |M*} fc=1 100Q of the estimator M n 
as described previously. We summarize our results in Table 5. We report that 
the obtained values v n are larger when 77 = 0.9 than when 77 = 1 . It is due to 
the noise effect. 

From now on, we fix the value of v n as in Table 5. With the same protocol as 
above and for all considered cases, we compute 1 000 other independent values 
{Mifj, , , „„„ of the estimator M n , for n = 50 000. Hence, we evaluate the 

l n J k=l,...,l 000 

empirical power of our testing procedure and the empirical first type error. 
Tables 6 and 7 provide the empirical results obtained by our test procedure in 
the experiments we deal with. We see that our testing procedure provides very 
good results even in the framework B, we obtain powers II equal to 1 since 
iV = 13 when we are in presence of noise. Otherwise, for N = 14 and rj = 0.9 
the powers of our testing procedure is a little bit degraded, but still remarkably 
good, since II = 0.7160 for a = 1% and II = 0.9440 for a = 5% in the framework 
B when p is the coherent-^/& state, the optimal N when we are in presence of 
noise is N = 13 with powers of test equal to 1. 

Table 5 

Empirical values of v n for r the vacuum state (case A) and r the coherent-3 state (case B) 







N = 15 
v = l 


N = 14 
n = 0.9 


N = 13 
77 = 0.9 


Case A 


u n (l%) 


0.0096 


0.4605 


0.2692 






0.0079 


0.2739 


0.1666 


Case B 


fn(l%) 


0.0028 


0.1453 


0.0798 




U n (S%) 


0.0022 


0.1100 


0.0713 
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Table 6 

Empirical values of the first type error a and the power of the test II over 1 000 runs for v, 
given in Table 5 for r the vacuum state ( case A ) 



a) p vacuum b) p single ph. c) p Schrodinger (7.-3 

an n 



N = I5,n = 1 


fn(l%) 


0.0150 


1.0000 


1.0000 




fn(5%) 


0.0570 


1.0000 


1.0000 


N = 14, n = 0.9 


fn(l%) 


0.0070 


1.0000 


1.0000 




f«(5%) 


0.0550 


1.0000 


1.0000 


N = 13, n = 0.9 


fn(l%) 


0.0130 


1.0000 


1.0000 




Vn(5%) 


0.0680 


1.0000 


1.0000 








Table 7 




Empirical values of the first type error a 


and the power of the test n over 1 000 runs for v n 




given in 


Table 5 for r the coherent-3 state (case B) 






a) p coherent-3 b) p coherent-\J% 


c) p Schrodinger C.-3 






a 


n 


n 


N = 15, 77 = 1 


Un(l%) 


0.0090 


1.0000 


1.0000 






0.0530 


1.0000 


1.0000 


N = 14, 77 = 0.9 




0.0120 


0.7160 


1.0000 




fn(5%) 


0.0680 


0.9440 


1.0000 


N = 13, n = 0.9 


fn(l%) 


0.0230 


1.0000 


1.0000 




fn(5%) 


0.0620 


1.0000 


1.0000 



5. Proof 

In this section, we give the proofs of the Theorem 1 and 2 derived in Section 3.4 
and the proofs of the Proposition 1 and 2 established in Section 3.3. 

5.1. Proof of the Theorems 

In the paragraphs below, we establish the results of the Theorem 1 and 2 derived 
in Section 3.4. We begin by Theorem 1-(1). 

Proof of Theorem 1-(1). Take r €]0,2], 77 = 1, the bandwidth Ni defined 
in the equation (18) and (p„ defined in (19), we have by Proposition 2-2 and 
Proposition 2-1 

V T (M„) < ^^iV 1 17/3 < Cy^l (24) 

V p (M n ) < ^l N ^ < Qvtpt (25) 

where Cy and Cy T are positive constants depending only on B, r and L. 

Under p = r defined in H , from Proposition 1, equations (24) and (19), we 
bound from above the first type error as follows: 

PA «, - 1] = P, [|M„, > a*) < « - « < -. 0, 



as n — * 00. 
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On the other hand, under p satisfying Hi(C,ip n ), we have the second type 
error bounded as follows: 

p p [n n = o] = p p [|m„| < c*t 2 n ] 



\M n - E p [M n }\ > \\p - r|| 2 - C*tl - B p [M n ] 



Moreover, under p satisfying Hi(C,ip n ), let C = C*(l+S). S e]0, 1[, Proposition 1, 
equations (29) and (23) imply that B p < xC*ip 2 for n large enough and then 



p p \n n = Q]<P P 



<Pr, 



\M n — E p [M n ] 
\M n -E p [M n ] 



< 4Cv < * 
" (<5C*) 2 " 2' 



for C* large enough. 



> 



> 



SC*/2 



C*tl-B p [M n 



□ 



Proof of Theorem l-(2). For r = 2, rj e]l/2,l[, the bandwidth N 2 defined 
in the equation (20) and ip^ defined in (23), we know by Proposition 2-2 and 
Proposition 2-1 that 



V T (M n ) < i^l!iV 2 - 4/3 e 32 ^ < C' v y ri 



(26) 
(27) 



where Cy and C v ^ are positive constants depending only on B, r, L and r\. 

On one hand, under p = r defined in Hq, by Proposition 1, equations (26) 
and (21), we have the first type error bounded as follows: 



P T [fl„ = l] = P T [\M n \ >C*t 2 n ] < 



E T [\M n \ 2 } 



V T [M n ] < ; Q 



C*Hi 



C* 2 



as n — > co. 

On the other hand, under p satisfying Hi(C,tp n ), the second type error is 
bounded as follows: 



P p [Q n = 0] = P p [\M n \<e*t 2 n ] 



< Po 



\M n - E p [M n \\ > \\p - rf - C*tl - Bp [M n ] 



From Proposition 1, equations (27) and (21), since we are under p satisfying 
Hi(C,<p n ), we get 
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M„ - E p \M n ]\ \\ P -t\\ 2 - C*tl - B p [M„ 



VVpWn 



\M n -E p [M n }\ > C-C*-C B 



c' 



fry 



(C-C*-C B ) 2 - 2' 



for C > C* and C* large enough. 



□ 



Proof of Theorem 2. In the case r G]0, 2[ and r\ G]l/2, 1[, for the bandwidth 
N$ solution of the equation (22) and for cp^ defined in (23), Proposition 2-2 and 
Proposition 2-1 give that 



;(M„) < ^Liv 3 - 4 /3 e 32 7 iV3 



„ 2(3r-14) „ 

<c v y n -N 3 3 <c v y n (28) 



V p (M n ) < — —N 3 ~ 



» »r-2/3 I67./V3 



(29) 



where Cy and C Vt are positive constants depending only on B, r, L and rj. 

Under p = r defined in Hq, from Proposition 1, equations (28) and (23), it 
follows that the first type error satisfies 

p tP . . i] = p, [|m„| > «a < « . M < _ 0, 

as n — ► 00. 

Furthermore, under p satisfying Hi(C, ip n ), the second type error is such that 



P p [n n = 0]=P p |M n | <C*£ 



|M„ - E p [M n }\ > \\p - r|| 2 - C*4 - i? p [M n ] 



From Proposition 1, equations (29) and (23), since we are under p satisfying 
Hi(C, ip n ), we deduce that 



p p [n n = o]<p p 



\M n -E p [M n ]\ > C-C*-C B 



VV P (M„) 

as n — > 00, for C > C* and C* large enough 



f7 AT 6 

°V JV 3 



< 



(C-C* - C B ) 2 



□ 



5.2. Proof of the Propositions 



In the following paragraphs, we give the proofs of Propositions 1 and 2 estab- 
lished in Section 3.3. 
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B p (M n ) = \E p [M n ] - \\p-r\\l\ = 



e p [M n ] ~ E \P3* ~ T i> h 

j,k>0 



It is easy to see that under p = r denned in Hq, we have B T (M n ) = 0. Under p 
satisfying Hi(C,tp n ), 

N-l 

B p {M n )= Y \Pj,k ~ Tj,k\ 2 ~ E \pj,k-Tj, k \ 
j,k=0 j,k>0 
oo j — 1 oo 

= 2 e e i^ fe ~ r ^ fc i 2 + e ~ r wi 2 

j — N k=0 j=N 

OO — 1 OO 

< 4 E E (i*.*i a + M 2 ) + 2 E 0*«i a + i^i 2 ) 

j=iV k=0 j=N 

Since r and p belong to the class 7Z(B, r, L) defined in (If), it implies 

oo j—1 



B p (M n ) <8L 2 J2J2 e- 2BU+k) + 4L 2 £ 



j=N k=0 



j=N 



POO POO 

< 8L 2 / (j - f )e- 2Bur/2 du + AL 2 / e - aB W"du 

JN JN 

< C B N 2 - r / 2 e- 2BNr/2 + C' B N x - r l 2 e- 2B{ - 2N ^'\ 

where Cb and C' B denote positive constants depending only on B, r and L. As 
2 r ' 2 > 1 for all r > 0, 

B p {M n ) < C B N 2 - r / 2 e- 2BNr/2 (l + o(l)), 



as TV ■ 



□ 



Proof of Proposition 2. By centering variables, we have M n — E p [M n ] 
F\ + F2 where 



1 N-l n , 



pV { — $ 



N-l 



FV.i — ,$ 



j,k=Q 1=1 



F\ ( — ,$ 



x E n 
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vV 



where IZe(z) stands for the real part of the complex number z. 

Then, V p (M n ) := E„[\M n - E p [|M„|]| 2 ] = £ p [|Fi| 2 ] + S P [|F 2 | 2 ] . First deal 
with the first term of the previous sum £^[1^1 1 



E„ 



\Fi\' 



n(n — 1) p 



N-l 

E 



j,k=0 

- E 



F v 



Y 



21. 



By the Cauchy-Schwarz inequality on the sum and as E[\X — E[X]\ 2 ] < E[\X\ 2 ] 

2 



E„ 



1 



N-l 



j,fc=0 



F\ [ — ,* 



A direct consequence of Remark 2 is 



N -A/i e Wr,N for r/G] 1/2,1 [, 



< 



(£1^^17/3 f Qr7? = l. 

By noticing |7£e(.z)| < \z\, the second term of the sum E p [Ffl is such that 



< -E„ 



x 



JV-l 



j,k=0 



y 



By the Cauchy-Schwarz inequality on the sum and as — i?[X]| 2 ] < E^Xl 2 ]: 

'jv-i 



E a 



W2Y 



< 



E*p 

J,k=0 



JV-l 

E 

U,k=0 



E 



Y 



<I> 



T 3,k 
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Figure 3. For r the vacuum state (case A): a) p the vacuum state, b) p the single photon 
state and c) p the Schrodinger cat-3 state for rj = 1, N = 15. 




Single photon 



Schrodinger cat 



Figure 4. For t the vacuum state (case A): a) p the vacuum state, b) p the single photon 
state and c) p the Schrodinger cat-3 state for n = 0.9, JV = 14. 



By Remark 2 and as \\pW2 , IMI2 < 1 we obtain under p satisfying H\(C, ip n ) 



8C1 



JV -2/3 e 16 7 n forjjG]l/2,l[, 



I -Fa 



< 



N 17 / 6 for n=l. 
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Figure 5. For r the vacuum state (case A): a) p the vacuum state, b) p the single photon 
state and c) p the Schrodinger cat-3 state for n = 0.9, N = 13. 
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Figure 6. For r the coherent-^ state (case B): a) p the coherent-3 state, b) p the coherent-^/6 
state and c) p the Schrodinger cat-3 state for n = 1, N = 15. 



If, N is such that the upper bound of i? p [|F2| 2 ] tends to with n, then this 
is the dominant term in the upper bound of V p (M n ). In addition, notice that 
under p = r defined in H" , we have E T [\F 2 \ 2 } =0. □ 
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Figure 7. Forr the coherent-3 state (case B): a) p the coherent-3 state, b) p the coherent-\^6 
state and c) p the Schrodinger cat-3 state for i) = 0.9, N = 14. 
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Figure 8. Forr the coherent-3 state (case B): a) p the coherent-3 state, b) p the coherent-^/6 
state and c) p the Schrodinger cat-3 state for n = 0.9, N = 13. 



6. Appendix 

The boxplots of the 1000 values {M„} ft _ 1 10Q0 of the estimator M n (for 
n = 50 000) in the case A and B are represented by Figures 3 to 5 and by 
Figures 6 to 8 respectively. 
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